% 隔壁 chemlab.m 的示例
% 改变相应的参数吧！

% 1. 化合反应 A+B->C
N = 3;
N_eq = 1;

% 初始物质的量，nA=1, nB=1, nC=0;
n(1,:)=[1,1,0];

Eq(1).nu =[-1,-1,1];
Eq(1).kp = 5;
Eq(1).kr = 1;

% 2. 分解反应
% 其余配置同化合反应，只需改变初始时各物质的的量
% 由于化学反应平衡，C将会自发分解

% 初始物质的量，nA=0, nB=0, nC=1;
n(1,:)=[0,0,1];

% 3. 中间体理论
% 这个例子展示了一个包括反应中间体的例子
% 该反应由两步组成, A+B->C, C->D, C是反应中间体
% C很难被生成，也很快会分解。尽管C的含量一直很少，但是C是反应的必要一环

N = 4;
N_eq = 2;

n(1,:)=[1,1,0,0];

% A+B->C
Eq(1).nu =[-1,-1,1,0];
Eq(1).kp = 1;
Eq(1).kr = 10;

% C->D
Eq(2).nu =[0,0,-1,1];
Eq(2).kp = 100;
Eq(2).kr = 1;

% 4. 催化
% 无催化路径：A+B->C
% 加入催化剂 D 后：
% A+D->E
% B+E->C+D
% 这两个路径比原路径更快
% 尽管催化剂D的量很少，也不被“消耗”，也不改变化学反应平衡，但是D提高了反应速率
% 这个例子受W站“催化”词条启发

N = 5;
N_eq = 3;

n(1,:)=[1,2,0,0.01,0]; %改变D的浓度以调整有无催化剂。

% A+B->C
Eq(1).nu =[-1,-1,1,0,0];
Eq(1).kp = 0.06;
Eq(1).kr = 0.01;

% A+D->E
Eq(2).nu =[-1,0,0,-1,1];
Eq(2).kp = 200;
Eq(2).kr = 100;

% B+E->C+D
Eq(3).nu =[0,-1,1,1,-1];
Eq(3).kp = 300;
Eq(3).kr = 100;

